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The phenomenon of drag reduction by polymer additives had been studied in simulations on 
the basis of non-Newtonian fluid mechanical models that take into account the field of polymer 
extension (conformation tensor) and its interaction with the velocity field. Drag reduction was 
found in both homogeneous and wall bounded turbulent flows. In the latter case it was shown 
recently that the notion of scale-dependent effective viscosity allows quantitative predictions of 
the characteristics of drag reduction in close correspondence with experiments. In this paper we 
demonstrate that also drag reduction in homogeneous turbulence is usefully discussed in terms of 
a scale-dependent viscosity. In other words, the essence of the phenomena under study can be 
recaptured by an "equivalent" equation of motion for the velocity field alone, with a judiciously 
chosen scale-dependent effective viscosity that succinctly summarizes the important aspects of the 
interaction between the polymer conformation tensor field and the velocity field. We will also clarify 
here the differences between drag reduction in homogeneous and wall bounded flows. 



I. INTRODUCTION 

The addition of long chained polymers to turbulent 
flows can result in a significant reduction in the drag 
Q I^JM The phenomenon had been discovered in 
1949 Q and had since attracted large amount of atten- 
tion, with much of the experimental literature reviewed 
and systematized by Virk Q; the amount of drag de- 
pends on the characteristics of the polymer and its con- 
centration, but cannot exceed an asymptote known as 
the "Maximum Drag Reduction" curve which is indepen- 
dent of the polymer's concentration or its characteristics. 
The understanding of this phenomenon had seen signifi- 
cant progress recently. A first step in forming a new un- 
derstanding were direct numerical simulations of model 
equations of viscoelastic flow^ both in wall bounded and 
in homogeneous turbulence ja, 0, ■ The Oldroyd-B and 
the FENE-P models first, and then simplified models like 
shell models and Burger's like models all ex- 

hibited drag reduction as a result of including the in- 
teraction between the velocity field and a second field 
representing the polymer (end-to-end) conformation ten- 
sor, see Figs. 1 and 2. In homogeneous turbulence drag 
reduction is exhibited as the increase in the root-mean- 
square (rms) velocity fluctuations at scales larger than 
the Lumley scale defined as the scale for which the eddy 
turnover time is of the order of the polymer relaxation 
time. The intermediate scale rms energy fiuctuations are 
suppressed due to transfer of energy to the polymers. In 
wall bounded turbulence drag reduction entails an in- 
crease in the mean velocity for a given pressure head, 
see Fig. ^ Here the Reynolds stress at the intermedi- 
ate scales is suppressed [ij]; we will argue however that 
there is a difference between the increase in the rms ve- 
locity fluctuations at large scales in homogeneous flows 
and the increase in mean velocity in wall bounded flows; 
the former disappears when the system size goes to in- 
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FIG. 1; Mean velocity profile of the FENE-P (dashed 
line) and of the Navier-Stokes equations (solid line) in wall 
bounded channel flow as a function of the reduced distance 
from the wall. The relative increase of the mean velocity (in- 
dicated by the asymptotic straight lines) is the phenomenon 
of drag reduction in wall bounded flows. 



flnity (for a flxed Lumley scale). In the latter case an 
increase in the mean velocity near the wall (small and in- 
termediate scales) does not disappear with increasing the 
system's size. This difference is fundamental to the dif- 
ferent symmetries at play, the Galilean invariance in the 
case of the wall bounded flow vs. translational invariance 
in the case of homogeneous flows. Nevertheless we will 
argue below that the two cases can be discussed in simi- 
lar physical terms. In a recent paper it was shown that 
drag reduction in wall bounded flows can be conveniently 
discussed in terms of a 'scale-dependent' effective viscos- 
ity. The aim of the present paper is to demonstrate that 
this notion is also useful in the context of homogeneous 
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FIG. 2: Energy spectrum of the SabraP model (line) and the 
Sabra model (dashed line with symbols) for v = 10~^. The 
relative increase of the energy spectrum at small values of n is 
the phenomenon of drag reduction in homogeneous turbulence 
in general and in shell models in particular, see Sect. Illll for 
details. 



turbulence. In doing so we aim at simplifying the theo- 
retical description, eliminating the explicit presence of a 
second field in the equations of motion, leaving the ve- 
locity field alone. The eliminated field, which represents 
the conformation tensor of the polymers, remains only as 
an effective viscosity in the equation of motion. Needless 
to say, this effective viscosity cannot be a number, since 
the amount of energy transferred from the velocity field 
to the polymer is strongly scale dependent; in homoge- 
neous turbulence this transfer achieves a maximum near 
the Lumley scale. In wall bounded flows the degree of 
interaction between the polymers and the velocity field 
is a strong function of the distance from the wall, and 
so is therefore the effective viscosity. Of course, in a full 
theory a scale-dependent scalar viscosity is not sufficient 
either, due to the anisotropy of the polymer end-to-end 
extension tensor. We would like to demonstrate however 
that at least in the model equations a surprising propor- 
tion of the essential physics can be captured in terms of 
a simple notion of a scale-dependent viscosity which sur- 
rogates the existence of the second field. This thinking 
goes back to some observations a few years ago regard- 
ing the importance of space-dependent viscosity even in 
the stability of laminar flows [l4lll5l|. In Sect. 2 we re- 
view the two-fleld models in which drag reduction had 
been demonstrated in numerical simulations. In Sect. 3 
we present the reduction to velocity- alone models with 
scale-dependent viscosity. In Sect. 4 we present a dis- 
cussion of the large system size limit and underline the 
difference between homogeneous and wall bounded flows. 
Sect. 5 is dedicated to a short summary and conclusions. 



II. SHELL MODEL FOR DRAG REDUCTION 
IN HOMOGENEOUS TURBULENCE 

Viscoelastic flows are represented well by hydrody- 
namic equations in which the effect of the polymer en- 
ters in the form of a "conformation tensor" i?(r, t) which 
stems from the ensemble average of the dyadic product of 
the end-to-end distance of the polymer chains J[,9, ,2QJ . A 
successful model that had been employed frequently in 
numerical simulations of turbulent channel flows is the 
FENE-P model. Flexibility and finite extendability of 
the polymer chains are reflected by the relaxation time 
r and the Peterlin function P{r,t) which appear in the 
equation of motion for R: 
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In these equations and p\ refer to the maximal and 
the equilibrium values of the trace R~f^. Since in most 
applications Pm ^ Po the Peterlin function can also be 
written approximately as P{r,t) « (1/(1 — aR^^) where 
a = p~^. In its turn the conformation tensor appears in 
the equations for the fluid velocity u{r^t) as an additional 
stress tensor: 
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vs^^u + V • r + F , (3) 
(4) 



Here Vs is the viscosity of the neat fluid, F is the forcing 
and Vp is a viscosity parameter which is related to the 
concentration of the polymer, i.e. Vp/ ~ $ where $ 
is the volume fraction of the polymer. We note however 
that the tensor fleld can be rescaled to get rid of the pa- 
rameter in the Peterlin function, R^p — aRafj with 
the only consequence of rescaling the parameter po ac- 
cordingly. Thus the actual value of the concentration is 
open to calibration against the experimental data. These 
equations were simulated on the computer in a channel or 
pipe geometry, reproducing faithfully the characteristics 
of drag reduction in experiments. It should be pointed 
out however that even for present day computers simu- 
lating these equations is quite tasking. It makes sense 
therefore to try to model these equations further. For 
the purpose of studying drag reduction in homogeneous 
systems one can derive a shell model whose simplicity 
and transparency are assets for analysis and simulations 
alike. In developing a simple model one is led by the 
following ideas. First, it should be pointed out that all 
the nonlinear terms involving the tensor field R{r, t) can 
be reproduced by writing an equation of motion for a 
vector field B{r,t), and interpreting i?^^ as the dyadic 
product BaBfj. The relaxation terms with the Peterlin 
function are not automatically reproduced this way, and 
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one needs to add them by hand. Second, we should keep 
in mind that the above equations exhibit a generahzed 
energy which is the sum of the fluid kinetic energy and 
the polymer free energy. Led by these consideration the 
following shell model was proposed in |^ ^ : 
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In these equations u„ and Bn stand for the Fourier am- 
plitudes u{kn) and B{kn) of the two respective vector 
fields, but as usual in shell model we take n = 0, 1, 2, . . . 
and the wavevectors are limited to the set A:„ = 2". The 
nonlinear interaction terms take the explicit form 



<^n{u,B) 



+ (1 + b)kn_2Un-2Bn-i 



(6) 



with b a parameter and the obvious extension to <i>„(u, u), 
^n{B, u) and $„(-B, B). In accordance with the general- 
ized energy of the FENE-P model 19, 20], also our shell 
model has the total energy 
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FIG. 3: The average value of the Peterlin function P{B) as 
a function of c computed in the SabraP model. The dashed 
line corresponds to P = 1. 



III. SCALE DEPENDENT EFFECTIVE 
VISCOSITY IN HOMOGENEOUS DRAG 
REDUCTION 



E = 




(7) 



The second term in the generalized energy contributes 
to the dissipation a positive definite term of the form 
iiyp/T^)P^iB)J2n \Bn\^- With lyp = the first of Eqs. 
(O reduces to the well-studied Sabra model of Newtonian 
turbulence. We therefore refer the model with i/p ^ as 
the SabraP model. As in the FENE-P case we consider 
c = Vp/vs to be proportional to the concentration of poly- 
mers. In it was shown that this shell model exhibits 
drag reduction, and the mechanism for the phenomenon 
was elucidated. Furthermore, it was shown in [Tlj that 
for large enough concentration, the Peterlin function can 
be disregarded (i.e. P « 1) and, consequently, the dy- 
namics of the system becomes concentration indepen- 
dent, i.e. we reach the MDR asymptote. This behavior 
of the Peterlin function is shown in Fig. |2| Following 
the above finding, we consider below the limiting case in 
which the concentration is large enough for the Peterlin 
function to be close to unity, P « 1. Finally, all the 
numerical simulations reported in this paper have been 
performed by using b = —0.2, i/g = 10^^ and a constant 
energy input given by: 
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for n = 1,2 and P„ = for n > 2. 



(8) 



Drag reduction in homogeneous turbulence is exhibited 
by a relative increase in the rms fiuctuations of the energy 
at large scales. We thus focus naturally on the energy 
spectrum e(fc„) = (u„w*). In the context of the shell 
model the phenomenon is demonstrated in Fig. |5] where 
e{kn) is shown for the given values of the parameters. 
The spectra for the pure Sabra model (line with symbols) 
and the coupled model (line) are compared for the same 
amount of power input per unit time. The discussion of 
the spectra revolves around the typical Lumley scale kc 
which is determined by the condition 



= 1/2 



{kc)k(. 



{P{B)) 



(9) 



For kn ^ kc the decay time r becomes irrelevant for the 
dynamics of P„. The nonlinear interaction between u„ 
and Bn at these scales results in both of them having 
the same spectral exponent which is also the same as 
that of the pure Sabra model. The amplitude of the 
Un spectrum is however smaller than that of the pure 
Sabra in the coupled case, since the P„ field adds to the 
dissipation. On the other hand, for kn ^ kc the P„ field 
is exponentially suppressed by its decay due to r, and the 
spectral exponent of u„ is again as in the pure Sabra. 
Drag reduction comes about due to the interactions at 
length scales of the order of kc which force a strong tilt in 
the Un spectrum there, causing it to cross the pure Sabra 
spectrum, leading to an increase in the amplitude of the 
energy containing scales. This is why the kinetic energy 
is increasing for the same amount of power input, and 
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FIG. 4: The spectrum of the energy dissipation k^e{kn) for 
the Sabra (solid Une with symbol) and the SabraP models 
(dashed line). Both models show the same maxima for about 
n ~ 15 which corresponds to the peak in the energy dissipa- 
tion. 
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FIG. 5: Time behavior of Sp, as defined in Eq. (IIOII . which 
represents the whole energy exchange from the Un field to Bn . 
Negative values of Sp means that energy is taken from u„. 
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hence drag reduction. In Fig. 2|we show the spectrum of 
energy dissipation k^e{kn)- This figure indicates that as 
far as the dissipative scale is concerned, it is not changed 
much by the coupling of the velocity field to the polymer 
field; both models show a maximum at n ~ 14 which is 
the dissipative scale. We now address the question how to 
recapture the same phenomenon in a model involving the 
velocity field alone but with a scale-dependent effective 
viscosity. We first reiterate that the field u„ loses energy 
in favor of the field Bn- Using Eq. ((Sj we can measure 
the energy transfer from m„ to Bn using the quantity: 



VpP{B) 
3t 



Re{iE„<$„(B,S)} . (10) 



This function measures the exchange between the kinetic 
energy S„u*u„ and the "polymer" or "elastic" energy 
S„i?*i3„. In Fig. El we show a snapshot of the depen- 
dence of the function Sp on time. The point to notice is 
that Sp is negative for most of the time. The Bn field 
drains energy from the velocity field, and we therefore can 
hope to be able to capture its role by an effective viscos- 
ity. Note however that the dynamics of Sp is strongly 
intermittent; this feature is common to the shell model 
and the full FENE-P model as observed in the DNS of 
the latter. We cannot hope to capture all the temporal 
complexity with the notion of effective viscosity, since 
the latter is an average notion. Nevertheless the essen- 
tial features will be shown to be reproduced. We will 
try to capture the effect of Sp in terms of an effective 
viscosity as follows: using (..) for the (time) average, we 
introduce the scale dependent effective viscosity Vf,{kn) 
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FIG. 6: The values of the eddy viscosity Ve(kn) defined in Eq. 
mil for P{B) = 1. Note that this quantity rises rapidly in 
the vicinity of the Lumley scale. 



Ve(kn) 



j{kn)) 



fc^e(fc„) 



(11) 



The quantity i^e(^n) is shown in Fig. its maximum is 
reached at n 6—7, a wavenumber which is not yet in the 
dissipative range. It is important to stress that z^e(fc„) is 
obtained by averaging over a complex and intermittent 
dynamical behavior of the viscoelastic shell model. It 
is therefore not obvious that the main characteristics of 
drag reduction can be obtained by simply replacing the 
viscoelastic terms $„(i?,i?) by a scale dependent effec- 
tive viscosity. We demonstrate that this is possible by 
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FIG. 7: The energy spectrum of the SabraP model (sohd hne) 
as compared with the energy spectrum of the Sabra model 
with the effective viscosity and a = 0.3 (symbols) and the 
Sabra model without effective viscosity (dashed line). 



using now the Sabra model with an extra viscous term 
given by Ve{kn)k'^Un- The new viscous term replaces, on 
the average, the effect of viscoelastic terms proportional 
to ^n{B. B). The equations of motion read: 

'liT ~ l^eikn)klun - I'sk'^Un + Fn ■ (12) 

We do not expect that Veikn) in the dynamics of the 
Sabra model, as proposed in Eq. (|12|l . will be exactly the 
object measured on the average defined in Eq. l(TT|l . We 
clearly must keep the functional dependence of i^e on , 
but we can allow a factor of proportionally that will take 
care of the difference between the dynamical intermittent 
behavior and the average behavior. We will therefore use 
the form a^'e(fc„), where a is a constant that can be op- 
timized to achieve a close correspondence between the 
two-field model and the effective one-field model. For 
a = we recapture the original Sabra model without ef- 
fective viscosity. We simulated the Sabra model with the 
effective viscosity [Eq. (|12ll ] for different values of a in the 
range (0, 1). Drag reduction was found in all cases. For 
a = 0.3 the energy spectrum turns out to be very close 
to the original SabraP model with the viscoelastic terms. 
In Fig. 13 we show the energy spectrum of the SabraP 
model and the energy spectrum of the Sabra model with 
effective viscosity for a = 0.3. In order to check that 
the result shown in Fig. [T] are due to a scale dependent 
dissipation, we have defined a scale independent viscosity 
V* as: 

V* — ("^p) 1^23) 

5jnfc^e(fc„) 



The definition of v* is similar to that given in Eq. 1)11(1 , 
i.e. V* is defined such that, by adding a viscous term 
to the Sabra model, the system on the average 
is losing the same amount of energy as in the case of 
viscoelastic flows. It turns out that in our case v* ~ 
2.5 X 10^^. By using this value for v* we have numerically 
integrated the Sabra model by adding a new viscosity 
equal to z/*, namely: 

Hr ^ 3*"("' ") ~ ^ Ki^n - l^sknUn + Fn (14) 

The corresponding energy spectrum is shown in Fig. |S1 
together with the energy spectrum for the Sabra model 
(viscosity and the Sabra model with the effective vis- 
cosity 0.3ve{kn)- As one can clearly see, an increase of 
the dissipation for all scales does not result in a drag 
reduction. Finally, we have computed the energy flux 
of the Sabra model with effective dissipation and com- 
pared it against the energy flux of the SabraP model. 
This comparison is exhibited in Fig. |3 where the solid 
line corresponds to the SabraP model and the symbols 
correspond to the Sabra model with effective dissipation 
0.3i^e(fcn). The two energy fluxes are equal in the inertial 
range up to wavenumber n ^ 7. 

The results illustrated so far support the conclusion 
that a scale dependent effective viscosity is able to re- 
produce most of the dynamics of viscoelastic terms and, 
in particular, the phenomenon of drag reduction. Let us 
remark once more, that it is the scale dependence of the 
effective viscosity which is able to properly reproduce the 
drag reduction. It is worthwhile to explain the mecha- 
nism of the action of the scale dependent viscosity, to 
understand its similarity to the action of the polymers. 
For fixed energy input, as in our case, drag reduction 
is shown as an increase of the rms fluctuations at scales 
larger than the Lumley scale. The scale dependent ef- 
fective viscosity increases the viscous terms fc^w„ in a 
particular range of scales, say for Uc < n < n2, where 
Uc = log2(fcc). The energy flux n„ in the system is given 
by the third order correlation function 

n„ ~ (u* u„_|_i). As shown in the Fig. |51 we 
can safely assert that the energy flux docs not change for 
n < Uc- The increase of viscosity at n — Uc produces a 
decrease of the energy at scale Uc- Thus, we expect to 
decrease with respect to the value observed in the Newto- 
nian case. Since n„ is not affected by the increase of the 
viscosity at n = Uc, we must conclude that the quantity 
Un-iUn should increase while decreases. This is the 
origin of the tilt in the An increase of u„ spectrum in the 
vicinity of n^. From a physical point of view, this picture 
is not different from the one discussed in where a simi- 
lar explanation for the drag reduction was invoked. Note 
that all that we need for the phenomenon to occur is that 
the increase in viscosity should start at the right scale. 
This scale is equivalent of the Lumley scale whose role in 
the viscoelastic case had been already emphasized. 

Finally, we discuss the effect of changing the concentra- 
tion on the effective viscosity. When {P{B)) > 1 the ef- 
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FIG. 8: Energy spectrum of the Sabra model for Us (line), 
the Sabra model with increased viscosity v* (symbols) and for 
the Sabra model with an effective viscosity 0.3i^e(A;n) (dashed 
line). 
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FIG. 9; Energy flux computed for the SabraP model (solid 
line) and the Sabra model with effective viscosity (symbols). 



fective viscosity depends on the Peterlin function, which 
in turn depends on the concentration c and on the re- 
laxation time r, cf. Eq. Figure [TUI displays the 
effective viscosity as a function of fc„ for four values of 
the concentration, c = 10~^, 10~^, 10 and 100. As the 
concentration decreases, the effective viscosity decreases, 
and its peak migrates to higher values of This mi- 
gration is simply due to the change in the Lumley scale, 
cf.. Eq. @. The decrease in the effective viscosity is due 
to the increase in {P{B)) shown in Fig. |2| Needless to 
say, these changes in the effective viscosity decrease the 
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FIG. 10; Effective viscosity for varying the concentration: 
c — 10~^ (circles), c = 10~^ (squares), c = 10 (triangles) and 
c = 100 (line). 



effect of drag reduction, as seen in experiments and sim- 
ulations: only large concentrations agree with the MDR 
asymptote. 



IV. THE LIMIT OF LARGE SYSTEM SIZE 

In this section we want to discuss the limit fco ^ 
while keeping fixed the scale and the shape of the effective 
viscosity. In other words, we study fco — > for fixed 
value of the Lumley scale kc- Note that we take kc much 
smaller than the dissipative scale and we keep constant 
the rate of energy input e. 

The discussion simplifies by considering the other typ- 
ical scale in our system, which is the Taylor microscale 
At, 



(15) 



In j9| it was shown that the conditions are optimal for 
drag reduction in our shell model when a dimensionless 
parameter fi = Xxkc, is of the order of unity. On the 
other hand drag reduction is lost when /i 1 or /x <C 1. 
Obviously, when fco ^ the overall kinetic energy in- 
creases as /cq ^^"^ while the denominator in Eq. (|15(l re- 
mains unchanged, being dominated by the viscous scale. 
Thus fco ^ leads to Xt oo, and we expect to lose 
drag reduction in that limit (for a fixed value of kc). This 
conclusion is supported by the results shown in Fig. 1111 
where we plot the ratio between the kinetic energy with 
the effective viscosity and the Newtonian kinetic energy 
for L = k^^ ^ oo. The case L = 1 corresponds to the 
previous sections. Note that for L large enough, the sys- 
tem exhibits drag enhancement. Physically, for very large 
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FIG. 11: Ratio of the kinetic energy for the Sabra model with 
scale dependent viscosity and the kinetic energy of ths Sabra 
model with fixed kinematic viscosity, for different values of 
L = k^^ . Note that for drag reduction to take place the ratio 
must be larger than 1. The position the maximum in the scale 
dependent viscosity is kept fixed while L — > oo. 



values of k^/ko the effective dissipation is just increasing 
the overall viscosity in the system and, therefore, no drag 
reduction can be observed. For drag reduction to occur 
we must have the Lumley scale close to energy contain- 
ing scales. Note, however, that "close" in our case means 
kc ^ 50 — 100 larger than the integral scale ko. 

It is interesting to compare our findings, which pertain 
to homogeneous systems, to drag reduction in turbulent 
boundary layers. The elastic layer in such flows (between 
the viscous layer and the Newtonian plug) has the pe- 
culiar distinction that y, the distance from the wall, be- 
comes the only important scale in the problem. It is both 
the energy containing scale and the Lumley scale at the 
same time. The former is clear; at distance y from the 
wall the most energetic eddies are of size y. The latter 
needs a bit of theory, and this is provided in . The 
upshot of the analysis there is that in the elastic layer the 
kinetic energy K{y) scales like K{y) ^ y^ jr^ . Thus the 
Lumley scale is also y. Accordingly, the phenomenon of 
drag reduction is totally indifferent to the physical size 
of the channel (or pipe). As long as the conditions for 
drag reduction hold at distance y from the wall, drag 
reduction will occur and will have a persistent effect on 
the mean flow independently of the outer scale. Even- 
tually, when y is large enough, K(y) may stop growing 
like the Lumley scale decreases, and we observe cross 
over back to the Newtonian log layer, albeit shifted to a 
larger value of a mean velocity profile. 

In summary, drag reduction phenomena in homoge- 
neous and wall bounded flows have a lot in common even 
though the effect disappears in the former when the sys- 
tem size goes to infinity. The essential physics is the 
proximity of the Lumley scale to the energy containing 
scales, which allows an effective interaction between the 



V. CONCLUSIONS 

The work presented in this paper supports two con- 
clusions. First, we demonstrated that drag reduction by 
polymers can be represented in terms of an effective scale 
dependent viscosity. One can use a theory in which two 
fields are explicitly presented, i.e. the velocity field and 
the polymer field. Then the viscosity remains Newto- 
nian, and the polymer conformation tensor acts as the 
additional sink of energy at the intermediate scales which 
are larger than the viscous scales but smaller than the 
Lumley scale. We can construct however effective mod- 
els in which only the velocity field is present, and replace 
the polymer field by an effective viscosity. This effec- 
tive viscosity will be different from the Newtonian one at 
the crucial scales at which the polymers are active, i.e. 
scales larger than the dissipative scales but smaller than 
the Lumley scale. With a properly chosen effective vis- 
cosity we can reproduce the results of the two-field theory 
qualitatively and even semi-quantitatively. Having done 
so, we reach a unified discussion of drag reduction by 
polymers in homogeneous and wall bounded flows. It is 
worth pointing out however that the unified discussion is 
deeper than the device of unified viscosity. Superficially 
drag reduction in homogeneous and wall bounded tur- 
bulence appear very different. In the former there is no 
mean fiow and drag reduction appears as an increase of 
the rms fluctuations of the large scales. In the latter drag 
reduction means the increase of the mean fiow velocity. 
Nevertheless in essence the phenomenon of drag reduc- 
tion in homogeneous and wall bounded fiows is basically 
the same: the polymers act to reduce the gradients at the 
intermediate scales. They partly laminarize the flow at 
the intermediate scales, and this allows the largest scales 
to attain higher rms fluctuation levels (in homogeneous 
flows) or higher mean velocity (in wall bounded fiows). 
To understand this further recall that for laminar flows 
the drag is a strongly decaying function of Re. Once tur- 
bulence sets in, the dramatic increase in eddy viscosity 
contributes to a drag which is much larger than the one 
that would obtain in a hypothetical laminar flow with 
the same value of Re. The addition of polymers allows 
one to bring the drag closer to the hypothetical laminar 
low value, and this is done by reducing the turbulence 
level at intermediate scales. Whether one prefers to de- 
scribe the quantitative aspects of this phenomenon using 
explicitly the polymer field or by employing an effective 
viscosity depends to a large extent on one's goals. We ex- 
pect that the concept of effective viscosity will be found 
equally useful in discussing drag reduction in other situa- 
tions, for example when microbubbles are used instead of 
polymers. The quantitative aspects of such a description 
need however to be worked out case by case, and this is 
our program for the near future. 
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